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Abstract. We consider the problem of finding a function defined on (0, oo) from a countable set of values of 
its Laplace transform. The problem is severely ill-posed. We shall use the expansion of the function in a series of 
Laguerre polynomials to convert the problem in an analytic interpolation problem. Then, using the coefficients of 
Lagrange polynomials we shall construct a stable approximation solution. Error estimate is given. Numerical results 
are produced. 
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1. Introduction 

Let L^(0, oo) be the space of real Lebesgue measurable functions defined on (0, oo) such that 

poo 

||/||£ s = / |/(x)| 2 e-^<oo. 
" Jo 

This is a Hilbert space corresponding to the inner product 

POO 

<f,9>= / f(x)g(x)e~ x dx. 
Jo 

We consider the problem of recovering a function / S L 2 p {0, oo) satisfying the equations 

poo 

Cf{p 3 ) = / e-»*f(x)dx = n j (DIL) 
Jo 

where Pj € (0, oo), j = 1, 2, 3, ... 

Generally, we have the classical problem of finding a function f(x) from its given image g(p) 
satisfying 

poo 

Cf{p)= / e-**f{x)dx = g(p), (1) 
Jo 

where p is in a subset to of the complex plane. We note that Cf(p) is usually an analytic function 
on a half plane {Re p > a} for an appropriate real number a. Frequently, the image of a Laplace 
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transform is known only on a subset u> of the right half plane {Re p > a}. Depending on the set u, 
we shall have appropriate methods to construct the function / from the values in the set 

{£f(p) : pecu}. 

Hence, there are no universal methods of inversion of the Laplace transform. 

If the data g(p) is given as a function on a line (— ioo + a, +ioo + a) (i.e., oj = {p : p = a + iy, y e 
R}) on the complex plane then we can use the Bromwich inversion formula ([26], p. 67) to find the 
function f{x). 

If u> C {p G R : p > 0} then we have the problem of real inverse Laplace transform. The right 
hand side is known only on (0,oo) or a subset of (0,oo). In this case, the use of the Bromwich 
formula is therefore not feasible. The literature on the subject is impressed in both theoretical and 
computational aspects (see, e.g. [2,3,10,16,18,22]). In fact, if the data g{p) is given exactly then, by 
the analyticity of g, we have many inversion formulas (see, e.g., [3,7,8,20,21,23]). In [3], the author 
approximate the function / by 

N 

/(t) - ^(^(eW))/^ 

fc=0 

where bk(a) are calculated and tabulated regularization coefficients and g is the given Laplace trans- 
form of /. Another method is developped by Saitoh and his group ([4,5,20,21]), where the function 
/ is approximated by integrals having the form 

poo 

wjv(i) = / g(s)e- st P N (st)ds, N=l,2,... 
Jo 

where P^ is known (see [5]). Using the Saitoh formula, we can get directly error estimates. 

However, in the case of unexact data, we have a severely trouble by the ill-poscdness of the 
problem. In fact, a solution corresponding to the unexact data do not exist if the data is nonsmooth, 
and in the case of existence, these do not depend continuously on the given data (that are represented 
by the right hand side of the equalities). Hence, a regularization method is in order. In [7], the 
authors used the Tikhonov method to regularize the problem. In fact, in this method, we can 
approximate uq by functions up satisfying 

/3up + £*£up=£*g, (3 > 0. 

Since C is self-adjoint (cf. [7]), the latter equation can be written as 

I3u + / ^ds = / e- st g( S )ds. 
Jo s + t J 

The latter problem is well-posed. 

Although the inverse Laplace transform has a rich literature, the papers devoted to the problem 
with discrete data are scarce. In fact, from the analyticity of £f(p),if£f(p) is known on a countable 
subset of oj C {Re p > a} accumulating at a point then £f(p) is known on the whole {Re p > a}. 
Hence, generally, a set of discrete data is enough for constructing an approximation function of /. 
It is a moment problem. In [15], the authors presented some theorems on the stabilization of the 
inverse Laplace transform. The Laplace image is measured at N points to within some error e. This 
is achieved by proving parallel stabilization results for a related Hausdorff moment problem. For a 
construction of an approximate solution of (DIL), we note that the sequence of functions (e~ PjX ) is 
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(algebraically) linear independent and moreover the vector space generated by the latter sequence is 
dense in L 2 (0, oo). The method of truncated expansion as presented in ([6], Section 2.1) is applicable 
and we refer the reader to this reference for full details. In [11, 13], the authors convert (DIL) into a 
moment problem of finding a function in L 2 (0, 1) and, then, they use Muntz polynomials to construct 
an approximation for /. 

Now, in the present paper, we shall convert (DIL) to an analytic interpolation problem on the 
Hardy space of the unit disc. After that, we shall use Laguerre polynomials and coefficients of 
Lagrange polynomials to construct the function /. An approximation corresponding to the non 
exact data and error estimate will be given. 

The remainder of the paper divided into two sections. In Section 2, we convert our problem into 
an interpolation one and give a uniqueness result. In Section 3, we shall give two regularization 
results in the cases of exact data and non exact data. Numerical comparisons with exact solution 
are given in the last section. 

2. A uniqueness result 

In this paper we shall use Laguerre polynomials 

X Jll 

L »w = 

We note that {L n } is a sequence of orthonormal polynomials on L^O, oo). We note that (see [1], 
[9], page 67) 



exp 

Hence, if we have the expansion 



then 



It follows that 



— J (1 - a)" 1 = £ £„(*)*"■ 

oo 

f(x) = a n L n{x) 
n-0 

/■oo / \ 00 

/ /(*) exp — - )(l-z)- 1 e- x dx = J2a n z n . 



f°° ( x \ 
y^a n z n = i f(x) exp I ) (1 - z y 1 dx. 

Put $f(z) = Yl™=o a nZ n , OLj = 1 - one has 

i.e., we have an interpolation problem of finding an analytic function <j>/ in the Hardy space H 2 (U). 
Here, we denote by U the unit disc of the complex plane and by H 2 (U) the Hardy space. In fact, we 
recall that H 2 (U) is the space of all functions <p analytic in U and if, (f> <G H 2 (U) has the expansion 

4>{z) = EfcLo a k zk then 

mWo = £m 2 = ^ / W 6 )\ 2 de. 
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We can verify directly that the linear operator $ is an isometry from L 2 onto H 2 (U). In fact, we 
have 

Lemma 1 Let f G L 2 (0, oo). Then Cf{z) is analytic on {z G C| Rez > 1/2}. If we have an 
expansion 



f — O-nLn 
n=0 

then one has $/ G H 2 (U) and 

oo 

ii*/iiW) = Ei a «i 2 = n/!i i?(o,oo)- 

n=0 

Moreover, If we have in addition that \[xf G Lp i/ien 

oo 

E™KI 2 <ll^/'lli ? . 



n=0 



Proof 

Putting F z (t) = e~ zt f(t), we have F z G L 2 (0,oo) for every Rez > 1/2. Hence £/(z) = 
Jo°° Fz{t)dt is analytic for i?ez > 1/2. From the definitions of Lp(0,oo) and H 2 (U), we have the 
isometry equality. Now we prove the second inequalities. We first consider the case /', /" in the 
space 

B = {g Lebesgue measurable on (0, oo)| \[xg G L 2 (0, oo)}. 
We have the expansion 

oo 

/ = O-nLn 
n=0 

where a n =< f,L n >. 

The function y — L n satisfies the following equation (see [17]) 



which gives 
It follows that 



na„ 



xy" + (1 - x)y' + ny = 
(xe~ x y')' + nye~ x = 0. 

/•OO 

/ f(x)nL n (x)e~ x dx 
Jo 

/•OO 

- / f(x)(xe- x L' n (x))'dx 
Jo 

/•oo 

/ f l (x)xe- x L' n {x)dx 
Jo 

{f(x)xe- x )'L n (x)dx 



f 

Jo 



(xf'(x) + f'(x) - xf'(x))L n (x)e- x dx 
= -<xf" +f -xf',L n > . 
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Since L n is an orthonormal basis, we have the Fourier expansion 

oo 

xf"+f'-xf' = J2(-™n)L n . 

n=0 

Using the Parseval equality we have 

oo 

< Xf" + f - Xf', f > = ^(-Mn)fln. 

n=0 

It can be rewritten as 



/ {xe- x f{x))'f{x)dx = -Y,nal. 

JO „ = n 



n=0 

Integrating by parts, we get 



/>oo - 

/ xe-*\f'(x)\ 2 dx = J2na 2 n . 

Ja n=0 



Now, for /' e B we choose (/&) such that f' k , /" k & B for every k — 1,2, ... and \fxf k (resp./fe) 
y/xf (resp./) in L 2 p as k — > oo. Assume that 

oo 
n=0 

Then we have 

„oc oo 
/ — -x| f'c v ,\|2 j \ v 

t fcn- 



/>oc °° 
JO r,-n 



n=0 

The latter equality involves for every N 

AT 



Y,™ln< l|Vi/fclli ? (O l0 o) ( 2 ) 



n=0 



Since fk — » / in L 2 as fc — » oo we have that fc — > oo, for each n. On the other hand, 

we have \/xf' k — > v 7 ^/' in L 2 as fc oo. . Therefore, letting k — > oo in (2) we get 



n=0 

Letting N — > oo in the latter inequality, we get the desired inequality. 
Using Lemma 1, one has a uniqueness result 

Theorem 1. Let pj > 1/2 for every j = 1,2,.... If 

XT 1 \ ^ 2 Pj - 1 

> h > — = oo 

, Pi , Pi 



then Problem (DIL) has at most one solution in L p (0, oo). 
Proof 

Let /i,/2 £ Z/p(0,oo) be two solutions of (DIL). Putting g = f\ — f 2 then g e ij;(0,oo) and 
Cg{pj) = 0. It follows that $3(1 — = 0, j = 1, 2, ... It follows that a 3 ■ — 1 — are zeros of 

<&g. We have <E>g e H 2 (U) and 

X>-kd=E^ + E ^ = °°- 

j=l Pj>l J 1/2< P] <1 Fd 

Hence we get $3 = (see, e.g. ,[19], page 308). It follows that g = 0. This completes the proof of 
Theorem 1. ■ 

3. Regularization and error estimates 

In the section, we assume that (pj) is a bounded sequence, pj 7^ pk for every j ^ k. Without 
loss of generality, we shall assume that p = 1 is an accumulation point of pj . In fact, if pj has an 
accumulation point po > 1 then, by putting f(x) = e~^ 0_1 ^ x /(a;) and p'j = pj — po + 1, we can 
transform the problem to the one of finding / e L 2 p (0, 00) such that 



f 

Jo 



e p i x f(x)dx = Hj, j = 1,2,... 



in which p'j has the accumulation point p = 1. In fact, in Theorem 2 below, we shall assume that 



1 



< a for every j — 1,2,..., where a is a given number. 



We denote by £ k m ^ (v) the coefficient of z k in the expansion of the Lagrange polynomial L m (y) (y 



(vi, v. m )) of degree (at most) to — 1 satisfying 

L m {v)(z k ) = v k , 1 < k < to, 
where z k = a k . If (f> is an analytic function on U, we also denote 

L m (4>) = L m {4>{z{), 4>{z m )). 

We define 

^MW - E 4 m) M* fc - 

0<fe<S(m-l) 

The polynomial L 6 m {v) is called a truncated Lagrange polynomial (see also [25]). For every g <E 
Lp(0, 00), we put 

T n g = (pi£g(pi),...,p n £g(pn)), 
Tg = (p 1 Cg(p 1 ),...,p n Cg( Pn ),...) ee°°. 

Here, we recall that a n — 1 — l/p n . We shall approximate the function / by 

F m = *- 1 L e m (T m f)= E A m \T m f)L k . 

O<fc<0(m-1) 
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We shall prove that F m is an approximation of /. Before stating and proving the main results, some 
remarks are in order. 

We first recall the concept of regularization. Let / be an exact solution of (DIL), we recall that 
a sequence of linear operator A n : £°° — ► L 2 (Q,oo) is a regularization sequence (or a regularizer) of 
Problem (DIL) if (A n ) satisfies two following conditions (see, e.g., [14], page 25) 

(Rl) For each n, A n is bounded, 

(R2) lim^oo || A, (T/) - /|| = 0. 

The number "n" is called the regularization parameter. As a consequence of (Rl), (R2), we can 

get 

(R3) For e > 0, there exists the functions n(e) and 5(e) such that lim e ^ n ( € ) = °°, lini e ^ 5(e) = 
and that 

\\A n(e) (p) - f\\ < 5(e) 

for every such that \\fi — TfW^ < e. 

In the present paper, the operator A n is <&~ 1 L m . The number e is the error between the exact 
data Tf and the measured data /i. For a given error e, there are infinitely many ways of choosing 
the regularization parameter n(e). In the present paper, we give an explicit form of n(e). 

Next, in our paper, we have the interpolation problem of reconstruction the analytic function 
<f> = $/ € H 2 (U) from a sequence of its values (<j)(a n ))- As known, the convergence of L m ((f>) to 
4> depends heavily on the properties of the points (a n ). The Kalmar- Walsh theorem (see, e.g. ,[12], 
page 65) shows that L m (4>) — > <j) for every <fi in C(U) for all <j> analytic in a neighborhood of U if and 
only if (a n ) is uniformly distributed in U, i.e., 

lim T /max|(z-ai)...(z-a m )| = 1. 

m^oo y \z\<l 

The Fejer points and the Fekete points are the sequences of points satisfying the latter condition (see 
[12], page 67). The Kalmar- Walsh fails if C(U) is replaced by H 2 (U) (see [25] for a counterexample). 
Hence, the Lagrange polynomial cannot use to reconstruct <p. In [12], we proved a theorem similar to 
the Kalmar- Walsh theorem for the case of H 2 (U). In fact, the Lagrange polynomials will convergence 
if we "cut off" some terms of the Lagrange polynomial. Especially, in [12] and the present paper, 
the points (a n ) are, in general, not uniformly distributed. 

In Theorem 2, we shall verify the condition (R2). More precisely, we have 

Theorem 2 Let a e (0, 1/3), let f G L 2 p (0, oo) and let p 3 > 1/2 for j = 1, 2, ... satisfy 

1-1 <a. 

Pj 

Put 0q be the unique solution of the equation (unknown x) 

2<r 1 - x 
l-o- ~ 

Then for 9 € (0, 9q), one has 

\\f-F m \\ 2 L 2 — >0 asm^oo. 
//, we assume in addition that \fxf <E L 2 (0, oo) then 

11/ - F m \\ 2 L2p < (1 + m6) 2 \\f\\ 2 L2p (^pj + ^HVi/'Hi^coo). 
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Proof 

We have in view of Lemma 1 



\\f-F m \\h p = £ i4 m) i 2 + E m 2 

O<fc<0(m-1) k>0(m-l) 



where S ™ = au — £ k (T m f ). We shall give an estimate for S ™ . In fact, we have 

m— 1 oo 

ii$/ - L m {T m f)r H2(U) = y, i4 m) i 2 + E m 



r } i 2 +E u - 12 

k— k—ni 

On the other hand, the Hermite representation (see, e.g. [12], page 59, [24]) gives 

w ro (*)($/)(C)dC 



$/(^)-L ro (T m /)(0) = -L / 
2™ Vet 



/ac/ w m (C)(C-^) 

where u m (z) = [z — a.\)...(z — a m ). Now, if we denote by = = ... = and 

4 m) = i 

4 m) = E <*h-<*jr (l<r<m), 

*/(C)dC 



Ps " 2m J du C s+1 



Wm(C) 

then we can write in view of the Hermite representation 

oo / k 



*f(z) - L m (T m f)(z) = E E(" 1 ) V --r« ^ 

fc=0 \r=0 / 

From the latter representation, one gets 

4 m) = E(-l) r ^U™>, 0<fc<m-l. 
Now, by direct computation, one has 



But one has 

k m ( e ie )| > (\e ie \ - - \a m \) > (1 - 

Hence 

< ||$/||^ (t /)(l-<7)- ro . 



We also have 



where Ct, = 



\®m— r 

r . Hence, we have 



(m) | < a m-r C r < jn-kyn 



m k\(m-k)\ 

/n i-e\ m 

i4 ro) i<(i+^)ii/iu ? (^J . 

From the latter inequality, one has in view of (3) 

\\f-F m \\i 2p <(i+ m enf\\iJf—) + j2 



k>m9 



For G (0, 6 ), one has 



Hence, we have 



0<- <- =1. 

1 — (7 1 — C7 



lim ||/-F m ||| 2 =0 



as desired, since on the one hand we have the comparison between an exponential with base b < 1 
and a power function and in the other hand the remain of a convergent series l a fc| 2 - 
Now if y/xf G L 2 p (0, 00) then one has since ^ > 1 and from Lemma 1 



00 1 00 1 



l\\2 
L 2 - 



m9 

k>m0 k=0 

This completes the proof of Theorem 2. ■ 

Now, we consider the case of non-exact data. In Theorem 3, we shall consider the condition (R3) 
of the definition of the regularization. Put 



Dr, 



max max 

l<n<m \\z\<R 



.(*) 



{z - a n )w' m (a n ) 
Let tp : [0, 00) — > R be an increasing function satisfying 

ip(m) > mD m , m=l,2, ... 

and 

m(e) = ^-^e" 3 / 4 )] - 1 

where [x] is the greatest integer < x. 

Theorem 3. Let a e (0, 1/3), let f,y/xf G ip(0,oo) and Zei pj > 1/2 /or j = 1,2, ... satisfy 

< a. 



1 



1 

Pj 

Put 80 be the unique solution of the equation (unknown x) 



l-<7 



1. 
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Let e > and let (fij) be a measured data of (£f(pj)) satisfying 

sup|i>j(£/(pj) - < e. 

j 

Then for 9 € (0, 6q), one has 

2m(e) 



11/ - ^Ll^Wl, < 2(1 + m{e)6f\\f\\l l + _^||^/'|||, + 2^. 



where z/| = pj/ij for j = 1, 2, ... 

Proof 

We note that 



(z-a 3 )<(«3)' 
It follows that 

m(y )\\oo < emD m . 

Hence 

\\L 6 jT m f) - L e m (S)\\ H 2 (u) < \\L m (T m f) - M*/)||oc < 
It follows by the isometry property of $ 

||/ - $- 1 L^(^)||| ? < 2||/ - F m ||| ? + 2\\^Ll(T m f) ^Ll^)f Ll 

\2|i/-||2 /2c 1 9 \ 2 _ . j 



< 2(1 + ^11/111,^— j 
+2e 2 ™ 2 ^. 

By choosing m = m(e) we get the desired result. ■ 

4. Numerical results 

We present some results of numerical comparison between the function f(x) given in L 2 (0, oo) 
and its approximated form F m as it is stated in Theorem 2. 

First consider the function f(x) = e~ x and its expansion in Laguerre series 



c ^ = E W (4) 

n>0 

So in the Hardy space H 2 (U), we have to interpolate the analytic function 

*/(*) = E = ^ (5) 



n>0 



by the Lagrange polynomial L m (T m f), interpolation defined by 



L m (T m f) [ 1 - - ) = Pi [ e-P>*e- x dx = (6) 
PiJ Jo Pi + 1 
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where — > 1 as i — > oo. 

On the interval (—1.8, +1.8) we have drawn in Fig.l the curves c~ x and its approximation 
L m (T m f )(x) for m = 10. If m = 12 there is divergence for our interpolation (Fig. 2) outside the 
interval (—1, +1). 



-1.5 -1 -0.5 0.5 



Fig. 1 

f(x)=exp(-x) 



approx. function 
exact function 



approx. function 
exact function 



-2 -1.5 -1 -0.5 0.5 1 1.5 



Fig. 2 

In our 2nd example we have chosen the function 

n>0 

In the Hardy space the function 



3^ V 3 I 3 + x 

n>0 y 



(7) 
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is approximated by the Lagrange polynomial L m (T m f) at the points (1 — } — 4p. )> K — > 1 as 
z — > oo. 

The Fig. 3 (resp.Fig.4) show the quite good convergence (resp. divergence) on the interval (—2.8, 2.8) 
with m = 4 (resp. m = 11). 

f(x)=exp(x/4) 

20 | 1 1 1 1 1 1 1 

approx. function 

exact function 




Fig.3 



approx. function 
exact function 



Fig. 4 

In both cases we have chosen 9q = 0.29 with a = 0.25 (6>o given by ^ _ — = 1, < a < ^). So 
in the 2nd case the truncated Lagrange polynomial is almost verified since 11 x 0.29 <~ 3.2. 
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